Última actualización: 27 mayo 2020
Las muestras de Fósforo y Nitrato se hace mediante campañas esporádicas con frecuencia de ~2 meses. En dichas campañas no se realizan aforos. A continuación se presentan los valores muestreados en el período 2004-2016 así como el caudal ajustado que de cada día (ver: tratamiento de caudales).
library(kableExtra)
qlt_41 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/rch41_calidad.RDS")
qlt_41$qobs = qlt_41$qobs*1.12
kable(qlt_41[,1:4], booktabs = TRUE,
caption = "Muestras de Fósforo y Nitrato en Paso Pache",
col.names = c("Fecha", "Nitrato (mg/l)","Fósforo (mg/l)", "Caudal (m3/s)")) %>%
kable_styling(c("striped", "bordered")) %>%
scroll_box(width = "600px", height = "200px")| Fecha | Nitrato (mg/l) | Fósforo (mg/l) | Caudal (m3/s) | |
|---|---|---|---|---|
| 1 | 2004-12-21 | 0.050 | 0.1600 | 5.83856 |
| 2 | 2005-01-25 | 0.230 | 0.1500 | 2.90640 |
| 3 | 2005-02-22 | 0.340 | 0.2100 | 6.28880 |
| 4 | 2005-04-19 | 0.150 | 0.1800 | 280.94304 |
| 5 | 2005-08-02 | 0.630 | 0.1200 | 37.36880 |
| 6 | 2005-11-15 | 0.050 | 0.4900 | 13.58112 |
| 7 | 2006-03-21 | 0.530 | 0.1500 | 38.51456 |
| 8 | 2006-07-04 | 0.320 | 0.0700 | 49.92736 |
| 9 | 2006-09-26 | 0.130 | 0.1100 | 23.18064 |
| 10 | 2007-01-23 | 0.510 | 0.0900 | 0.89824 |
| 11 | 2007-06-05 | 0.550 | 0.4400 | 15.77968 |
| 12 | 2007-09-11 | 0.570 | 0.1800 | 40.75680 |
| 13 | 2007-11-20 | 0.050 | 0.1100 | 10.61648 |
| 14 | 2008-01-22 | 0.850 | 0.2400 | 2.28032 |
| 15 | 2008-05-20 | 0.050 | 0.1200 | 3.38352 |
| 17 | 2009-03-17 | 0.790 | 0.2200 | 79.80336 |
| 18 | 2009-04-28 | 0.060 | 0.0800 | 1.36528 |
| 19 | 2009-08-26 | 0.400 | 0.1900 | 33.30208 |
| 20 | 2010-05-19 | 0.070 | 0.1600 | 4.41392 |
| 21 | 2010-07-20 | 0.220 | 0.3600 | 490.98000 |
| 22 | 2010-09-28 | 0.320 | 0.1400 | 24.90320 |
| 23 | 2010-11-16 | 0.470 | 0.1400 | 11.40832 |
| 24 | 2011-01-20 | 0.014 | 0.7030 | 1.97568 |
| 25 | 2011-03-17 | 0.090 | 0.2580 | 2.17056 |
| 26 | 2011-05-05 | 0.035 | 0.1390 | 4.50688 |
| 27 | 2011-07-21 | 0.330 | 0.2430 | 98.70336 |
| 28 | 2012-01-04 | 0.014 | 0.3090 | 2.21536 |
| 29 | 2012-03-08 | 0.130 | 0.2250 | 2.84144 |
| 30 | 2012-05-10 | 0.015 | 0.1580 | 2.71936 |
| 34 | 2013-01-10 | 0.230 | 0.2400 | 20.43440 |
| 35 | 2013-03-14 | 0.260 | 0.3200 | 19.37488 |
| 36 | 2013-07-25 | 0.430 | 0.2790 | 20.93280 |
| 37 | 2014-01-23 | 0.350 | 0.2220 | 17.04752 |
| 38 | 2014-03-20 | 0.290 | 0.2000 | 18.72416 |
| 39 | 2014-05-08 | 1.000 | 0.2070 | 22.21520 |
| 40 | 2014-07-10 | 0.360 | 0.1880 | 116.85856 |
| 41 | 2014-09-25 | 0.370 | 0.1920 | 24.79456 |
| 42 | 2014-11-13 | 0.020 | 0.0973 | 24.13824 |
| 43 | 2015-01-22 | 0.380 | 0.2000 | 4.70848 |
| 44 | 2015-03-19 | 0.016 | 0.2920 | 2.16160 |
| 46 | 2015-07-23 | 0.020 | 0.1290 | 2.17840 |
| 47 | 2015-09-24 | 0.200 | 0.1700 | 6.73680 |
| 48 | 2015-11-19 | 0.055 | 0.1450 | 4.25040 |
qlt_23 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/rch23_calidad.RDS")
qlt_23$qobs = qlt_23$qobs*1.65
kable(qlt_23[,1:4], booktabs = TRUE,
caption = "Muestras de Fósforo y Nitrato en Paso de los Troncos",
col.names = c("Fecha", "Nitrato (mg/l)","Fósforo (mg/l)", "Caudal (m3/s)")) %>%
kable_styling(c("striped", "bordered")) %>%
scroll_box(width = "600px", height = "200px")| Fecha | Nitrato (mg/l) | Fósforo (mg/l) | Caudal (m3/s) | |
|---|---|---|---|---|
| 1 | 2011-01-20 | 0.014 | 0.5260 | 0.11880 |
| 5 | 2012-01-04 | 0.014 | 0.1300 | 0.11880 |
| 6 | 2012-03-08 | 0.300 | 0.2090 | 6.53235 |
| 7 | 2012-05-10 | 0.080 | 0.0752 | 0.68310 |
| 8 | 2012-06-28 | 0.100 | 0.0905 | 2.29515 |
| 9 | 2012-09-20 | 0.340 | 0.0975 | 47.77905 |
| 10 | 2012-11-08 | 0.120 | 0.1140 | 8.29620 |
| 11 | 2013-01-10 | 0.190 | 0.0838 | 1.43550 |
| 12 | 2013-03-14 | 0.140 | 0.0745 | 1.85790 |
| 13 | 2013-07-25 | 0.310 | 0.0869 | 5.28000 |
| 14 | 2014-01-23 | 0.150 | 0.0892 | 6.35910 |
| 15 | 2014-03-20 | 0.110 | 0.1000 | 3.15975 |
| 16 | 2014-05-08 | 0.160 | 0.0748 | 5.35260 |
| 17 | 2014-07-10 | 0.450 | 0.0566 | 18.97335 |
| 19 | 2014-09-25 | 0.130 | 0.0906 | 7.84245 |
| 21 | 2014-11-13 | 0.180 | 0.1300 | 8.65260 |
| 22 | 2015-01-22 | 0.190 | 0.0780 | 1.82985 |
| 23 | 2015-03-19 | 0.020 | 0.0830 | 0.11880 |
| 25 | 2015-07-23 | 0.020 | 0.0450 | 0.41085 |
| 26 | 2015-09-24 | 0.110 | 0.0830 | 4.80975 |
| 27 | 2015-11-19 | 0.094 | 0.0500 | 0.87450 |
| 28 | 2016-01-21 | 0.130 | 0.1270 | 0.11880 |
| 29 | 2016-04-07 | 0.110 | 0.1010 | 2.46675 |
| 30 | 2016-06-30 | 0.090 | 0.0700 | 32.92080 |
| 31 | 2016-10-13 | 0.083 | 0.1570 | 2.32320 |
La siguiente figura presenta la relación entre el caudal unitario (lps/km2) y la producción unitaria de fósforo (mg/s/km2). Ha sido obtenida a partir de los muestreos de calidad de agua en Paso Pache y Paso de los Troncos. Los valores de caudal has sido tomados de las correspondientes estaciones de caudal.
library(plotly)
area_41 = 4896
area_36 = 2744
area_23 = 687
qlt_41$qobs = 1000*qlt_41$qobs/area_41
qlt_23$qobs = 1000*qlt_23$qobs/area_23
qlt_41$NO3_mgxs = qlt_41$NO3_mgxs/area_41
qlt_23$NO3_mgxs = qlt_23$NO3_mgxs/area_23
qlt_41$PT_mgxs = qlt_41$PT_mgxs/area_41
qlt_23$PT_mgxs = qlt_23$PT_mgxs/area_23
fit41_PT = lm(log(qlt_41$PT_mgxs)~log(qlt_41$qobs))
fit41_NO3 = lm(log(qlt_41$NO3_mgxs)~log(qlt_41$qobs))
fit23_PT = lm(log(qlt_23$PT_mgxs)~log(qlt_23$qobs))
fit23_NO3 = lm(log(qlt_23$NO3_mgxs)~log(qlt_23$qobs))
fit41_PT = round(c(exp(fit41_PT$coefficients[1]), fit41_PT$coefficients[2]),2)
fit41_NO3 = round(c(exp(fit41_NO3$coefficients[1]), fit41_NO3$coefficients[2]),2)
fit23_PT = round(c(exp(fit23_PT$coefficients[1]), fit23_PT$coefficients[2]),2)
fit23_NO3 = round(c(exp(fit23_NO3$coefficients[1]), fit23_NO3$coefficients[2]),2)
qlt_41$fit_PT = fit41_PT[1]*qlt_41$qobs^fit41_PT[2]
qlt_41$fit_NO3 = fit41_NO3[1]*qlt_41$qobs^fit41_NO3[2]
qlt_23$fit_PT = fit23_PT[1]*qlt_23$qobs^fit23_PT[2]
qlt_23$fit_NO3 = fit23_NO3[1]*qlt_23$qobs^fit23_NO3[2]PT <- plot_ly(data = qlt_41, x = ~qobs, y = ~PT_mgxs, type="scatter",
mode = 'markers',colors="blue", width=2, name="Paso Pache")
PT <- PT %>% add_lines(data=qlt_41,x = ~qobs, y = ~fit_PT,
type = 'scatter',
mode = 'line',
line = list(color = 'blue', width = 2),
name=paste0("PT = ",fit41_PT[1],"Q^",fit41_PT[2]))
PT <- PT %>% add_markers(data=qlt_23,x = ~qobs, y = ~PT_mgxs,
type = 'scatter',
mode = 'markers',
marker = list(color = 'red', width = 2),
name="Paso de los Troncos")
PT <- PT %>% add_lines(data=qlt_23,x = ~qobs, y = ~fit_PT,
type = 'scatter',
mode = 'line',
line = list(color = 'red', width = 2),
name=paste0("PT = ",fit23_PT[1],"Q^",fit23_PT[2]))
PT <- layout(PT, xaxis = list(title="Caudal (lps/km2)",range=c(-1,2),type = "log"),
yaxis = list(title="PT (mg/s/km2)", range=c(-2,1),type = "log"),
legend = list(x = 0.01, y = 1),
width=500,
height=500
)
PTLa siguiente figura presenta la relación entre el caudal unitario (lps/km2) y la producción unitaria de Nitrato (mg/s/km2). Ha sido obtenida a partir de los muestreos de calidad de agua en Paso Pache y Paso de los Troncos. Los valores de caudal has sido tomados de las correspondientes estaciones de caudal.
NO3 <- plot_ly(data = qlt_41, x = ~qobs, y = ~NO3_mgxs, type="scatter",
mode = 'markers',colors="blue", width=2, name="Paso Pache")
NO3 <- NO3 %>% add_lines(data=qlt_41,x = ~qobs, y = ~fit_NO3,
type = 'scatter',
mode = 'line',
line = list(color = 'blue', width = 2),
name=paste0("NO3 = ",fit41_NO3[1],"Q^",fit41_NO3[2]))
NO3 <- NO3 %>% add_markers(data=qlt_23,x = ~qobs, y = ~NO3_mgxs,
type = 'scatter',
mode = 'markers',
marker = list(color = 'red', width = 2),
name="Paso de los Troncos")
NO3 <- NO3 %>% add_lines(data=qlt_23,x = ~qobs, y = ~fit_NO3,
type = 'scatter',
mode = 'line',
line = list(color = 'red', width = 2),
name=paste0("NO3 = ",fit23_NO3[1],"Q^",fit23_NO3[2]))
NO3 <- layout(NO3, xaxis = list(title="Caudal (lps/km2)",range=c(-1,2),type = "log"),
yaxis = list(title="NO3 (mg/s/km2)", range=c(-2,1),type = "log"),
legend = list(x = 0.01, y = 1),
width=500,
height=500
)
NO3